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Recent review article by S. Samuel "On the speed of gravity and the Jupiter/Quasar 
measurement" published in the International Journal of Modern Physics D, 13 (2004) 
1753, provides the reader with a misleading "theory" of the relativistic time delay in 
general theory of relativity. Furthermore, it misquotes original publications by Kopeikin 
and Fomalont & Kopeikin related to the measurement of the speed of gravity by VLBI. 
We summarize the general relativistic principles of the Lorentz-invariant theory of prop- 
agation of light in time-dependent gravitational field, derive Lorentz-invariant expression 
for the relativistic time delay, and finally explain why Samuel's "theory" is conceptually 
incorrect and confuses the speed of gravity with the speed of light. 
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1. Introduction 

Exact mathematical solution of the light geodesic equation in time-dependent grav- 
itational field of arbitrary moving bodies has been constructed in a series of our 
publications 1 ' 2 - 3 . This solution predicts that a light particle (radio wave) is de- 
flected by the gravitational field of the moving body from its retarded, with respect 
to observer, position. The retarded position of the light-ray deflecting body origi- 
nates from the Lienard-Wiechert solution of the linearized Einstein equations and 
can be found by solving the retarded-time equation which is a null characteristic of 
the gravitational field. We proposed relativistic VLBI experiment to measure this 
effect of retardation of gravity by the field of moving Jupiter via observation of light 
bending from a quasar and successfully completed this experiment in September 
2002. The experiment testing observational phase-referencing technique with sev- 
eral reference calibrators is described in 5 and the results of the main experiment 
of September 2002 are published in Astrophysical Journal 6 . 

Samuel's paper 7 is an attempt to review the results of the experiment by mak- 
ing use of a linear Lorentz transformation of the static Shapiro time delay. This 
approach is physically insufficient for conceptual analysis of the experiment which 
requires matching of the first and second order effects in the relativistic theory of 
the time delay based on the Lienard-Wiechert solution of the Einstein equations. 
This is the main reason for the erroneous conslusions about the nature of the exper- 
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iment mispresented in 7 . In the present paper we outline the basic equations of the 
complete Lorentz-invariant theory of the time delay and show mistakes in Samuel's 
linearized "theory" 7 originating from his first publication 8 . 

2. Retardation of Gravity and the Lienard-Wiechert Potentials 

We denote the barycentric coordinates of the solar system as x a = (x°,a; J ), 
where x° = ct, x 1 — x, and c is the fundamental speed limit. The metric ten- 
sor g^ v = r]^ + hpv, where r}^ = diag(— 1, 1,1,1) is the Minkowski metric, and h a p 
describes gravitational field of the solar system in the linearized post-Minkowskian 
9 approximation^ We impose the harmonic gauge conditions d v h^ v — l/2d^h = 0, 
where d v = d/dx v denotes a partial derivative with respect to coordinate x" , and 
h = rfh^. 

Linearized Einstein equations in harmonic coordinates assume the following form 

.i* +v .)»r-^(r--i^T), (i) 

where V = d/dx l , is the stress-energy tensor of matter generating the gravita- 
tional field, and T = rf v Ty, v . In what follows, we take the stress-energy tensor in the 
standard form of moving point-like masses with Dirac's delta-functions as amtter's 
support 10 . In linear approximation the gravitational field of the solar system bodies 
moving along their orbits with acceleration can be found from Eq. (1) separately 
for each body. For Jupiter, the perturbation of the metric tensor, h^ u , is obtained 
as the Lienard-Wiechert solution of the Einstein equations (1) which yields 

= 2GMj2u»u" + Tr 

c 4 r a u a ' 1 ' 

Here Mj is the mass of Jupiter, u a = j(s)(c, vj(s j) is its four-velocity with vj(s) = 
dx(s)/ds, 7(s) = (1 — v 2 j{s)/c?)^ 1 / 2 and the null-cone distance, r a — x a — Xj(s), 
and Jupiter's worldlinc Xj(s) — (cs,xj(s)) are functions of the retarded time s, 
determined as a solution of null cone equation rj^r^r" — 0, that is 

s = t — —\x — xj(s)\ . (3) 
c 

Eq. (3) describes the propagation of Jupiter's gravity field from Jupiter to the point 
x with the fundamental speed c b . The letter c used here is the speed of gravity 
and should not be confused with the speed of light (the speed of a radio wave 
from a quasar measured by VLBI) as it has been mistakenly done by Asada 12 . 



a Greek indices run from to 3. Roman indices run from 1 to 3. The Greek indices are rised and 
lowered with the Minkowski metric. Bold letters denote spatial vectors. Repeated indices mean 
the Einstein summation rule. Euclidean dot and cross products of two vectors are denoted as a • b 
and a X b respectively. 

b Because Eq. (3) is the null cone of the gravitational field of moving Jupiter we used 6,11 a symbol 
c g in there instead of c to mark the parameter which we fit to observations in the data processing 
procedure. 
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This fundamental conclusion of general relativity can be tested experimentally by 
observing gravitational interaction of light from a background source (quasar, star) 
with the time-dependent field of moving Jupiter and/or other solar system bodies 

4 



3. Electromagnetic Phase 

Radio interfcromctry (VLBI) measures the phase ip of the wave front coming from 
a quasar. The phase is determined from the cikonal equation for electromagentic 
field (the quasar radio wave) in the geometric optics approximation of Maxwell's 
equations 10 

gTd^ip - , (4) 

where in the linear approximation g^ v = rf" v — . Solution of the eikonal equation 
is obtained by iterations after substitution of the Lienard-Wiechert gravitational 
potentials, Eq. (2), into Eq. (4) 16 



v 

if = (fia + - 
c 



k a x a + — ^ (k a u a )ln<P 



(5) 



where $ = —k a r a , ipo is a constant, v is the radio frequency, k a = (l,fe) is a null 
vector associated with the propagating radio wave, and the unit vector k is directed 
from the quasar to VLBI station c . 

Solution for eikonal given by Eq. (5) describes propagation of a plane elec- 
tromagentic wave scattered on the time-dependent gravitational potential of a 
point-like mass (Jupiter) moving with time-dependent velocity. We abandoned the 
acceleration-dependent terms in Eq. (5) because of their smallness. Influence of the 
acceleration-dependent terms has been analysed in 2 . We notice that Eq. (5) is 
Lorentz-invariant because it consists of the products of four-dimensional vectors re- 
maining invariant with respect to the Lorentz transformations of the both Maxwell 
and Einstein equations. In particular, the argument of the logarithmic fucntion is 
$ = ~k a r a = r — k ■ r, where r = \r\ = \pr ■ r, r = x — xj(s), and depend on 
the retarded position of Jupiter, xj(s), calculated at the retarded time s defined by 
Eq. (3) of the null characteristic of the gravitational field of Jupiter. Thus, precise 
interferometric measurement of phase ip of the electromagnetic wave scattered by 
the gravitational field of a moving gravitating body (Jupiter) allows to determine 
the null characteristic of the gravitational field and measure the ultimate speed of 
propagation of gravity which is expected to be equal to the speed of light in vacuum. 
This is because the measurable phase ip ~ ln$(s) is a logarithmic function of the 
retarded time s as evident from Eq. (5). 

It is important to notice that the null characteristic of the gravitational field 
connecting the point of observation, x, and the retarded position of Jupiter, xj(s), 



c Notice that vector k a does not coincide with the null characteristic of the gravitational field given 
by the interval x a — Xj(s) (see Fig. 1). 
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can not be confused with the null characteristic of the electromagnetic field prop- 
agating along the space-time direction defined by the null vector k a because this 
vector is not parallel to the null vector r a — x a — x a (s) (see Fig. 1). Unfortunately, 
the four-dimensional Minkowskian diagram of the experiment shown in Fig. 1, was 
not understood by Asada 12 and Samuel 8,7 who confused the null characteristic of 
the quasar radio wave and that of the gravitational field of Jupiter. 



4. The Lorentz-invariant Theory of Time Delay 

The Lorentz-invariant relativistic time delay equation, generalizing the static 
Shapiro delay, can be obtained outright from equation (5). We note that the phase 
ip of electromagnetic wave, emitted at the point Xq = (ci , x ) and received at the 
point x a = (ct,x), remains constant along the wave's path 10 . Indeed, if A is an 
affinc parameter along the path, one has for the phase's derivative 

dip _ dp_dx£_ 

dx dx- d\ ' ( ' 

due to the orthogonality of the light rays and their wave fronts. Eq. (6) means that 
ip(x a (Xj) = const, in accordance with our assertion. Equating two values of the 
phase at the event of emission, Xq , and that of observation, x a , and separating the 
time and space coordinates one obtains from (5) 

1, , . IGMj 1 - c~ x k ■ vj , , . 
t - t = -k ■ (x - x ) 3 In (r-k-r) , (7) 

where both the distance r = x — xj(s), r = \x — xj(s)\, and the retarded time 
s are defined by the gravity null-cone equation (3). The time delay (7) of light 
propagating through the gravitational field of an arbitrary moving body was derived 
first by Kopeikin and Schafer 2 who solved equations for light geodesies in the 
retarded Lienard-Wiechert gravitational field of the body. Klioncr 13 also obtained 
this expression for the case of a uniformly moving body by making use of the 
Lorentz transformation of the Shapiro delay (that is, both the Maxwell and Einstein 
equations) from static to moving frame. Samuel's paper 7 represents, in fact, rather 
convoluted linearized analogue of Klioncr's painstaking calculation 13 . 



5. Relationship Between Our Theory and Samuel's Calculations 

Let us approximate Eq. (5) by introducing two angles & and 9 between vector k and 
unit vectors p = R/ R and I = r/r correspondingly, where vectors R = x — xj(t) 
and r = x—xj(s) connect the point of observation x = x{t) with the present, xj(t), 
and retarded, x.j(s), positions of Jupiter, respectively. By definition cos 6 = k ■ p 
and cos# = k ■ I. The product $ = —k a r a = r(l — cos (9) ~ r6 2 /2 for small angles. 
Hence, the phase variation caused by space-time difference Sx a — (cSt, 5x) between 
two VLBI antennas is 

*-!(M^ + ^f), (8) 
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where 56 = n ■ Sx/r with n = I x (fe x I) as the impact vector of the light ray with 
respect to the retarded (due to the finite speed of gravity) position of Jupiter xj{s). 
Undetectable terms of order v,j/c have been neglected. The quantity St — t 2 — t\ is 
the measurable VLBI time delay and Sx = x 2 — x\ = B is a baseline between two 
VLBI stations. Since VLBI stations measure the same wave front, S(p — 0. Thus, 
for sufficiently small angle 9 Eq. (8) yileds 

St = c- 1 k-B + A , (9) 



c J r 9 

where we neglected small terms of the second order of magnitude. 

Eq. (10) can be also derived from Eq. (7) 11 . It exactly coincides with the leading 
order term in Eq. (3.9) of Samuel's paper 7 and allows to establish the following 
relationship between Samuel's notations for the Earth-Jupiter distance Rej and 
the "observable angle" 9 b s d , and our notations for the null-cone distance r and 
the angle 9. Specifically, neglecting the Earth finite-size effects in our calculations, 
one obtains Rej = r and 6 a b s = 9. In what follows, we prefer to keep on using our 
notations. 

Eq. (10) is the excess time delay caused by the scattering of the radio wave 
from the quasar on the gravitational Lienard-Wiechert potential of moving Jupiter. 
Though the terms depending on Jupiter's velocity do not show up in Eq. (10) ex- 
plicitly, they are surely incorporated in it implicitly through the retarded position 
of Jupiter, xj(s), entering both the distance r, the angle 9, and the impact param- 
eter vector n. Such implicit dependence of the relativistic expressions on velocity 
of gravitating body is typical for the, so-called, post-Minkowskian approximation 
scheme of general relativity operating with the retarded Lienard-Wiechert solutions 
of the Einstein equations 9 ' 14 . Thus, Eq. (10) describes dynamical situation since 
(because of the orbital motion of Jupiter) the retarded quantities, r, 9, and n, 
are not reduced to their static counterparts after their post-Newtonian expansion 
around the time of observation t. 

The post-Newtonian expansion of the retarded position of Jupiter, xj(s), in Eq. 
(10) around time of photon's arrival to observer, t, yields 

xj(s) = xj(t) + vj{t){s -t) + 0(s - t) 2 . (11) 

The difference s—t is calculated by solving the gravity null-cone Eq. (3). Substituting 
this solution to Eq. (11) yields 

xj(s) = xj(t) - -vj(t)R + 0(s - tf , (12) 

c 



d These quantities appear in Eq. (1.3) in Samuel's paper 7 along with the "distance of the closest 
approach" £ spontaneously without rigorous mathematical description so that their meaning is 
fuzzy. 
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where R = R\, and R = x — xj(t) is the vector lying on the hypersurface of 
constant time t and connecting the point of observation, x, and the present position 
of Jupiter, xj(t). 

The post-Newtonian expansion given in Eq. (12) originates from the retarded 
nature of the Lienard-Wiechert gravitational potentials and, thus, describes the 
effect of the retardation of gravity 4 . The presence of the retardation of gravity effect 
through the retarded position of Jupiter in the Lorentz-invariant Eq. (5) and its 
small-angle approximation Eq. (10) reflects the causal property of gravitational field 
which is a consequence of its finite speed. Causality and Lorentz-invariance of the 
gravitational field are tightly connected fundamental concepts and the measurement 
of the causal nature (retardation) of gravity in the Fomalont-Kopeikin experiment 
is equivalent to the proof that the gravitational field is Lorentz-invariant and vice 
versa. 

Samuel also uses the post-Newtonian expansion of the retarded coordinate of 
Jupiter in Eq. (5.3) of his paper 7 . He believes that the post-Newtonian expansion 
" arises because the position of Jupiter changes as the quasar signals travel from the 
Jupiter region to Earth". Jupiter does move as the quasar signal travels towards 
observer but the distance traveled by Jupiter is proportional not to the difference 
between the time t* of the closest approach of the quasar signal to Jupiter and 
the time of observation, that is t* — t, but to the difference between the retarded 
time s and the time of observation t, that is s — t, as clearly follows from Eq. (11). 
Therefore, the fuzzy concept of " the Jupiter region" , which is repeatedly used by 
Samuel without rigorous mathematical definition, is, in effect, the retarded position 
of Jupiter defined by the Lienard-Wiechert solution of the gravity field equations. 
This consideration makes it evident that the origin of the post- Newtonian expansion 
in Eqs. (5.3)-(5.6) of Samuel's paper 7 is caused by the speed of gravity which 
propagates from moving Jupiter towards observer as well as the quasar signal does. 
This point was emphasized in our papers 4,16 . Samuel overlooked the gravitational 
physics of the Jupiter-quasar experiment because of his approximate and, hence, 
insufficient solution of the problem of propagation of light rays in time-dependent 
gravitational fields. He was able to integrate the light-ray geodesies only for the case 
of the small-angle approximation « 1) when the time of the closest approach, 
t* , of the quasar signal to Jupiter is comparable with the retarded time s along the 
null characteristic of Jupiter's non-stationary gravitational field. It is for this reason 
that Samuel confused the time t* and the retarded time s and replaced the concept 
of the propagation of gravity from the retarded position of Jupiter by the concept 
of the propagation of light from "the Jupiter region" c . 

Substitution of the post-Newtonian expansion (12) to Eq. (10) yields 

A = A S + A R . (13) 



c Similar mistake has been done also by Will in 15 who used insufficiently elaborated c g - 
paramctrization of the Einstein gravity field equations. 
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Here 

A s = - (4GMj/c 3 Q)(N ■ B) (14) 

is the static Shapiro time delay caused by Jupiter's gravitational field at the time 
of observation, where N = p x (k x p) is the impact vector of the light ray with 
respect to the present position of Jupiter xj(t), and 



AGMj 
R c 4 i?9 2 



2{N ■ vj)N - (K x (vj x K) 



(15) 



is the post-Newtonian correction to the Shapiro time delay due to the the finite speed 
of gravity in the gravity null-cone equation (3). Eq. (15) is the same as Eq. (4) from 
6 . It describes the first post-Newtonian vj/c correction to A s and can be detected 
because of the amplifying factor ~ I/O 2 . Samuel's Eq. (5.6) is an approximate 
form of our Eq. (15). We notice that Samuel's notation for the angle 9i = 9 in 
our notations. Physical origin of the post-Newtonian correction A# to the static 
Shapiro time delay Ag is due to the Lorentz-invariant nature of the gravitational 
field caused by its finite speed of gravity as follows from Eq. (11). 

In his first paper 8 Samuel incorrectly assumed that the experiment directly 
compared the radio position of the quasar with the optical position of Jupiter, and 
that the direction of Jupiter was determined by " sunlight that has been reflected off 
of Jupiter" (see the second paragraph in section III of Samuel's paper 8 describing 
figure 1 which is similar with figure 2 of Samuel's paper 7 ). This assumption would 
correspond to direct measurement of the angle 9 and hence no vj/c terms would 
be observed since they are not evident in Eq. (10). This explains why Samuel has 
erroneously decided "that the v/c effects are too small to have been measured in the 
recent experiment involving Jupiter and quasar J0842+1845" 8 . The experiment, 
however, monitored the position of the quasar as a function of the atomic time by 
the arrival of the quasar's photons at the telescope, while the Jupiter's position 
entering the time delay Eq. (7) was determined separately by fitting VLBI data for 
the quasar to a precise JPL ephemeris, evaluated at the same atomic time as the 
arrival of a photon via standard transformations from ephemeris time to atomic time 
6 . The result of our fitting procedure was that Jupiter deflects light from its retarded 
position xj(s) but not from its present position xj(t). Thus, the difference A — As 
was measured and the vj/c correction A# was determined within precision of 20% 6 . 
The measurement of the post-Newtonian correction (15) to the Shapiro time delay 
(14) is direct demonstration that gravity does propagate with the same speed as the 
speed of light. Samuel's claim that "Fomalont and Kopeikin's announcement that 
the speed of gravity is the speed of light to within 20% has no content" is based on 
his inability to distinguish between gravitational and electromagnetic effects in the 
post-Newtonian expansion of the Lorentz-invariant time delay Eq. (7) predicting 
that any moving body deflects light by its gravitational field from the retarded 
position in accordance with the causal (Lorentz-invariant) nature of gravity. 



7, 2008 3:14 WSPC/INSTRUCTION FILE samuel-comm 



8 Sergei M. Kopeikin 

6. Discussion and Further Particular Comments 

6.1. On the Ideology of the Experiment 

The paper by Samuel 7 is an attempt to protect his misleading calculations 8 pub- 
lished in 8 which conceptual inconsistency was revealed in MM 6 . 17 . Unfortunately, it 
does not provide either new mathematical details or more deep insight to the prob- 
lem. By making use of a linearized Lorentz transform of a spherically-symmetric 
gravitational field Samuel succeeds in calculation of the first few terms of the differ- 
ential Shapiro time delay caused by moving Jupiter that approximates the original 
result by Kopeikin 4 (see Eq. (7)) given in terms of the retarded time, s, connect- 
ing the point of observation, x = x(t), and the retarded position of a light-ray 
deflecting body (Jupiter), xj(s). Our derivation of Eq. (7) makes it clear that the 
electromagnetic signal from a quasar is observed at the time, t, and Jupiter de- 
flects it at the retarded time, s = t — r/c, where r — \x(t) — xj(s)\ is the radial 
coordinate of Jupiter with respect to observer directed along the null characteristic 
of the gravitational field defined by the null cone equation r\ jlv r tl r v = 0. Both the 
retarded coordinate of Jupiter, xj(s), and the retarded time, s, originate from the 
Lienard-Wiechert solution of the linearized Einstein equations which is a hyperbolic 
(wave-type) D'Alembert equation. Lorentz-invariant theory of the propagation of 
light rays through the time-dependent field of the gravitational retarded potentials 
reveals that the light particle (photon) is deflected by moving Jupiter when it is lo- 
cated at the retarded position xj(s) due to the finite speed of gravity. Experimental 
testing whether Jupiter deflects light from its orbital position taken at the retarded 
time s due to the finite speed of gravity or at the time of observation t, is a direct 
probe of the numerical value of the speed of gravity which must be equal to the 
speed of light according to Einstein. This is the key idea of the experiment which 
has been put forward in our publication 4 and practically tested in September 2002 
6 . Unfortunately, the einsteinian gravitational physics of the experiment is greatly 
misunderstood and conceptually misrepresented both in Samuel's papers 8 ' 7 and in 

12,15,18 f 

6.2. What Was Observed and Tested in the Experiment 

Samuel prefers to re-express our original result for time delay (7) in terms of the 
"observable angle" 9 b s = 6 in terms of "the distance of the closest approach 
The vertex of this angle is at the point of observation, x(t), and it has two legs - one 
leg is directed in the sky towards the quasar and another one is directed towards 
retarded position of Jupiter, x.j(s). If both Jupiter and the quasar were observed 
simultaneously in radio or in optics, then, the legs composing the "observable angle" 



Critical discussion of the formally correct, but conceptually misleading point of view on the 
gravitational physics of the process of light scattering by the gravitational field of a moving body 
presented in 18 , requires more elaborated mathematical technique than that presented in the 
present paper, and will be given somewhere else. Some details are available in 19 . 
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9 would be formed by the null characteristics of electromagnetic signals emanating 
correspondingly from the quasar and from Jupiter to the observer. However, Jupiter 
was not observed by VLBI in the Fomalont-Kopeikin experiment at all, only the 
quasar was - so that only one leg of " the observable angle" is the null characteristic 
of electromagentic field. Jupiter affected the propagation of the quasar radio wave 
by its gravitational field and the effective position of Jupiter in the sky (the second 
leg of "the observable angle") was determined in the process of the data analysis of 
the residual phase of the quasar radio wave. The time delay Eq. (7) tells us that the 
direction in the sky connecting the observer and Jupiter is the null characteristic of 
the gravitational field of moving Jupiter. Thus, Samuel's "observable angle" 6 b s = 
is made of the null lines one of which belongs to the electromagnetic field (radio 
wave from the quasar) but another one belongs to the gravitational field of moving 
Jupiter, which deflects the radio signal of the quasar not instantaneously but with 
the retardation to comply with the Lorentz-invariant symmetry of the Einstein 
equations s . 

The goal of our measurement was to distinguish between the two hypothesis: 

(1) Einsteinian gravity. Jupiter deflects light by its gravitational field from its 
retarded position x.j(s) as shown in our time-delay Eq. (7); 

(2) Newtonian gravity. Jupiter deflects light by its gravitational field instan- 
taneously and its position in the time delay Eq. (7) must be taken at the 
time of observation, t. 

In effect, wc did not measure any angles in the sky directly during the time of 
observation as Samuel seems to believe. What was measured is the residual phase 
Sip of the quasar radio wave, allowing us to determine the difference between times 
of arrival of the same front of the quasar radio wave to two VLBI stations 6 . This 
difference depends on the position of Jupiter on its orbit and our goal was to prove 
that the observed differential time delay is affected by the gravitational field of 
Jupiter acting from its retarded position, xj{s), but not from its position, xj(t), 
taken at the time of observation. The difference between the two hypothesis gives 
rise to the post-Newtonian correction Eq. (15) to the Shapiro time delay (14) which 
had to be zero if the speed of gravity would be infinite. We proved experimentally 6 
that the first hypothesis is correct within 20% (that is the post-Newtonian correction 
A_r ^ 0) which means that the speed of gravity has the same value as the speed of 
light. 

Samuel misinterprets the physical origin of the observed retardation-of-gravity 
effect in the time delay of light because of his insufficiently elaborated mathematical 
solution of the problem of light propagation in the time-dependent gravitational 
field of a moving body. While our approach starts from the basic principles and 
solution of the gravity field equations, Samuel's development deals only with the 
final product of the theory - the time delay Eq. (9), which incorporates effects 



g See Fig. 1 for graphical illustration of this point. 
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of the propagation of both gravitational and electromagentic fields. When making 
Lorentz transformation of the time delay Eq. (9) Samuel has missed the point that 
this transformation transforms not only the electromagnetic but the gravitational 
field as well. Thus, the experiment effectively measures the discrepancy between the 
Lorentz-invariant symmetry of the gravitational field with respect to the Lorentz 
transformation symmetry of the electromagnetic field which could arise due to the 
possible difference between the speed of gravity and light. The experiment has 
confirmed the Lorentz invariance of the gravitational field because no difference 
between the speeds of gravity and light was found within 20% 6 . 

6.3. On the Relativistic Jargon and the Symbol c 

Einstein's theory of general relativity predicts that an electromagnetic wave is de- 
flected by the gravitational field of moving body (Jupiter) from its retarded position, 
xj(s), taken at the retarded instant of time s = t — r/c, where c is conceptually 
the fundamental speed of propagation of gravitational field because it enters the 
Lienard-Wiechert gravitational potentials. For historical reasons, it is customary to 
call the symbol c as "the speed of light". We emphasize however that this term is 
misleading when one considers the Lienard-Wiechert solution of the Einstein grav- 
ity field equations since in this case c is the speed of gravity characterizing the 
degree of opening of the null cone made up of the characteristics of the gravita- 
tional field. This interpetation of the symbol c entering the retarded coordinate of a 
moving gravitating body remains true in the case of the physical process of scatter- 
ing of electromagnetic wave on the Lienard-Wiechert gravitational potential of this 
body. This point was misinterpeted by Asada 12 , Will 15 , and Carlip 18 . Samuel 7,8 
seems to be also confused with the existing conventional terminology and is about 
to interpret any physical effect depending on c as associated with propagation of 
light because c is always "the speed of light" by his definition h . Einstein's gravity 
field equations contain "the speed of light" constant c but it does not mean that 
the measurement of this fundamental constant for gravitational field (as it is done 
in 6 ) is reduced (or equivalent) to the measurement of the physical speed of light 
even if the light is used for such measurement. Various fasets of the symbol c and 
the pitfalls related to the physical interpetations of various experiments have been 
discussed recently in n > 20 . 

6.4. On the Point of the Closest Approach 

Samuel's calculations are also based on the assumption that an electromagnetic 
signal is effectively deflected at the point of its closest approach to the moving 
Jupiter From the first glance this assumption looks plausible but turns out to 

h The same jargon is used by Asada 12 and Will 15 that prevents them to see the true nature of 
the retardation of gravity effect measured in the Fomalont-Kopcikin experiment 6 . 
'The same assumption was also used by Will 15 
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be generally incorrect in exact mathematical solution of the problem shown in the 
present paper as well as in our previous works 1 t 2 > 3 AA6,2i an( j the work of Klioner 
13 . First of all, gravitational field is long-ranged and the process of gravitational 
light deflection (electromagentic wave scattering) does not take place in one event. 
Relativistic time delay of light is the integral effect and, strictly speaking, it can 
be expressed in terms of the impact parameter of the light ray, taken at the time 
of the closest approach t* , only approximately (see 1:2,3 for more detail). But even 
in the case of the approximation the point of the closest approach of light to the 
moving body loses its physical content when the Lorentz transformation is applied 
to transform the time delay from the static to a moving frame. Advanced calculation 
(missed in 7,s ) makes it evident that the integral effect of the light-ray time delay 
is coming out from the retarded position of the light-ray deflecting body 2 ' 3,4 which 
is in accordance with the Lienard-Wiechert solution of the Einstein gravity field 
equations describing propagation of the gravitational field J . 

Recently Klioner 13 has calculated the deflection and the time delay of light in the 
gravitational field of a uniformly moving body by making use of the Lorentz trans- 
formation technique. In contrast to the linearized Samuel's "theory" 7,8 , Klioner's 
calculation takes into account all velocity-dependent terms in the Lorentz transfor- 
mation which significantly supersedes Samuel's consideration 7 ' 8 . Klioner's paper 
13 delivers an independent proof that Jupiter must deflect light from its retarded 
position Xj(s) in accordance with the Licnard-Wicchcrt solution of the gravity field 
equations 2,3,4,n . The effect of the retardation of gravity in the case of a uniformly 
moving body is equivalent to the property of the gravitational field to be Lorentz- 
invariant and to propagate with the fundamental speed c which, at the same time, 
is the parameter of the Lorentz transformation of the Einstein gravity field equa- 
tions. Our experiment has tested this property of gravity by observing that Jupiter's 
gravitational field deflects the quasar radio wave out of Jupiter's retarded position 
lying on the null characteristic of the gravitational field which connects the VLBI 
station and Jupiter and is given by the retarded time Eq. (3). 

Different mathematical technique 1,2,3,21 also demonstrates that the point of 
the closest approach of photon to the moving body plays no physical role because 
the time of the closest approach, t*, drops out of the final equation for the light 
bending and/or relativistic time delay. What matters is the retarded time s — t—r/c 
of the Lienard-Wiechert solution of Einstein's equations, where c is the speed of 
propagation of gravity from the moving body to the point of observation. The time 
of the closest approach, t*, approximates the retarded time s = t — r/c, if the 
impact parameter of the light ray to the moving body is small enough (see 1,4 
for more detail). This mathematical approximation is physically possible because 
the gravitational field propagates with the same speed as light in general theory 



'We emphasize once again that there is no light propagating from the retarded position of the 
light-ray deflecting body in the time delay Eq. (7) to observer as it was erroneously deduced by 
Asada 12 , it is gravity which propagates. 
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of relativity. In fact, Samuel 8 ' 7 unknowingly used this approximation to replace 
the retarded time s with the time of the closest approach t* which has led him to 
the confusion of the propagation of gravity effect with that of the propagation of 
the quasar signal. His statement that the "Shapiro time delay difference is due to 
relatively short- distance effects is correct but these short- distance effects occur near 
the retarded position of Jupiter xj(s) taken at the retarded time s while the time of 
the closest approach t* is an approximation which, if it is used improperly, simply 
misinterpets the gravitational physics of the Jupiter-quasar experiment. 

6.5. On the Minkowski Diagram of the Jupiter- Quasar Experiment 

Samuel's graphic analysis of the Fomalont-Kopcikin experiment is shown in figures 
1 and 2 of the paper 7 . Unfortunately they do not grasp the relativistic spirit of 
the Jupiter-quasar experiment. The Minkowski diagram is much more adequate for 
picking up the idea of the measurement of the speed of gravity with VLBI. Analysis 
of the experiment in terms of the Minkowski diagrams is given by Kopeikin in 
the paper 11 and in the proceedings of the 14th Midwest Relativity Meeting 19 . 
Minkowski diagram of the experiment shown in Fig. 1 of the present paper clearly 
demonstrates the origin of the observed retarded position of Jupiter as caused by the 
Lorentz- invariant nature of the gravitational field and its finite speed of propagation. 
A null direction connecting observer and Jupiter in the relativistic time delay Eq. 
(7) is that corresponding to the null characteristic of the gravitational field. This is 
because of two facts: (1) it originates from the retarded Lienard-Wiechert solution 
of the gravity field equations, and (2) there is no radio or light wave emitted or 
reflected by and propagating from Jupiter towards observer during the time of the 
experiment. Thus, the measurement of the direction of the null characteristic of 
the gravitational field through the retarded coordinate of Jupiter, xj(s), through 
the relativistic time delay is a direct confirmation of Einstein's prediction that 
gravitational field has the same speed of propagation as the speed of light. The 
reader should not interpret the retarded position of Jupiter shown in figure 2 of 
Samuel's paper 7 as caused by propagation of sunlight reflected from Jupiter as 
Samuel 8,7 erroneously used to believe. Unfortunately, this was either overlooked or 
ignored by referees of papers 8 ' 12 despite of our persistent attempts to draw their 
attention to this, inconsistent with observations, fact. As explained in 6,11 the radio 
emission of Jupiter or the sunlight reflected by Jupiter can not be, and was not, 
observed at any VLBI station due to the specific technical limitations of VLBI. 
Anyone who discusses the nature of the Jupiter-quasar experiment must learn how 
VLBI operates and detect radio signals before making any statements about the 
nature of the measured effect of the retardation of gravity. The retarded position 
of Jupiter, measured in the experiment through the best fit of the observed VLBI 
time delay to its theoretical value, is due to the finite speed of gravity and reflects 
the fundamental fact that the gravitational field is Lorentz-invariant and its null 
characteristics coincide with the null characteristics of the electromagnetic field 
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withing 20% 6 . 

6.6. On the speed of gravity parameter c g 

Einstein's theory of general relativity does not require introduction of any parameter 
for the speed of gravity c g because it is Lorentz invariant and the speed of gravity is 
equal to the speed of light. All general-relativistic effects associated with the speed of 
gravity can be easily and unambiguously identified through the retarded positions 
of the gravitating bodies in the post-Minkowskian solutions of the gravity field 
equations and equations of motion. Will 22 has decided to introduce the parameter 
c g for the speed of gravity to facilitate discussion of the relativistic effects associated 
with the speed of gravity from those caused by the " speed of light" c. Introduction 
of c g to the Einstein equations is not trivial and requires to retain all differential 
relationships of general relativity for any value of c g . This task was not fulfilled in 

22 

We have found 11 the c s -parametrization of Einstein's equations which preserves 
all differential and algebraic properties of those equations for any value of the param- 
eter c g . This c ff -parametrization requires introduction of a global unit vector field, 
V a , which goal is to keep the Lorentz invariance of the gravity field equations. The 
difference between c g and c leads to one kind of observable effects while existence 
of spatial components, V, of the vector field V a leads to apperance of the preferred 
frame effects (the PPN parameters a\, <x 2 ) which have been strongly limited by 
pulsar timing and lunar laser ranging observations 22 . Hence, in the paper 11 we did 
not analyzed the preferred frame effects caused by V and worked in the coordinate 
system, where V a = (1,0,0,0), to concentrate on the discussion of the retardation 
of gravity effect 6 . Therefore, spatial vector field V does not appear explicitly in our 
derivation of the relativistic time delay equation. However, correct transformation 
of our time delay equation, derived in 11 for the case of c g =^ c, from one frame 
to another requires accounting for the preferred frame vector filed V, which was 
not done by Samuel in deriving his equation (7.2) in 7 . Thus, Samuel's criticism of 
my paper 11 is completely unfounded and is based on his misunderstanding of the 
mathemtical formalism which has been worked out in 11 . 
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Fig. 1. Two null cones related to the experiment are shown. The gravity field cone describes the 
retarded Lineard-Wiechert solution of the Einstein equations and reflects the Lorentz invariant 
nature of the gravitational field. The light cone shows propagation of light from the quasar. The 
relativistic perturbation of a light ray measured by an observer takes place when the gravity cone of 
Jupiter passes through the observer. The VLBI experiment measures the Minkowski dot product 
<I> between two null vectors at the point of observation directed to the quasar and to Jupiter 
respectively. Had Jupiter not been detected at the retarded position on its world line the speed of 
gravity were not equal to the speed of light, and the Einstein principle of relativity for gravitational 
field would be violated. The experiment did not find any violation of the Lorentz-invariance of the 
gravitational field and confirmed that the speed of gravity equals to the speed of light within 20%. 



